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Abstract: The influence of large permanent dipoles on molecular orbital 
tomography via high-order harmonic generation (HHG) is investigated in 
this work. It is found that, owing to the modification of the angle-dependent 
ionization rate resulting from the Stark shift, the one-side-recollision 
condition for the tomographic imaging can not be satisfied even with the 
few-cycle driving pulses. To overcome this problem, we employ a tailored 
driving pulse by adding a weak low-frequency pulse to the few-cycle laser 
pulse to control the HHG process and the recollision of the continuum 
electrons are effectively restricted to only one side of the core. Then 
we carried out the orbital reconstruction in both the length and velocity 
forms. The results show that, the orbital structure can only be successfully 
reproduced by using the dipole matrix elements projected perpendicular to 
the permanent dipole in both forms. 
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1. Introduction 

The recent development in strong-field physics has opened wide perspectives for observing the 
structure and ultrafast electron dynamics in atoms and molecules with Angstrom and attosecond 
resolutions fl]45l. One of the most fascinating topics at the frontier of this area is the molecular 
orbital tomography (MOT) based on high-order harmonic generation (HHG), which was first 
proposed by J. Itatani et al in {6 1. By applying this method to a series of aligned molecules, they 
accomplished a tomographic reconstruction of the highest occupied molecular orbital (HOMO) 
of N2. Since then, the MOT has attracted a great deal of attention for its great significance in 
uncovering the basic properties of molecules and chemical reactions J7ffi3). In most of the 
works, the target molecular orbitals are symmetric with either gerade or ungerade symmetry, 
while the MOT of nonsymmetric orbitals is more difficult. It has been discussed by Van der 
Zwan et al [14| that the MOT of orbitals with arbitrary symmetry requires one to restrict the 



recollision of the continuum electrons to only one side of the core in the HHG process 1 15 
[16). Therefore, they suggested using extremely short few-cycle laser pulses to ensure the one- 
side-recollision and the molecular orbital of HeH 2+ was successfully reconstructed in their 
simulation. Besides, one can also consider using an co + 2(0 two-color laser pulse for the MOT 
of nonsymmetric orbitals (7][l7j . 

In addition to the arbitrary symmetry, another important characteristic of the nonsymmetric 
orbital is the permanent dipole. As a large permanent dipole would lead to an unnegligible Stark 
shift of the orbital energy, it would have influence on the MOT in two aspects. Firstly, it will 
cause additional laser-induce Stark phase acquired in the HHG of the nonsymmetric orbital. 
As the phases of field-free dipole matrix elements are needed in order to reconstruct molecular 
orbitals, one has to find a way to subtract the Stark phase from experimental measurements. 
The method is found in p8) , in which the authors derived an analytical expression of the first- 
order Stark phase only depending on the measurable quantities. Secondly, previous works have 
shown that the Stark shift of the orbital energy can efficiently modify the angular dependence of 
the ionization rate and the distribution of the ionized electrons [ 19-21 1. This effect will make it 
more complicated to control the continuum electron for one-side-recollision and to reconstruct 
the orbital. This influence has however not been discussed before to the best of our knowledge 
and requires more concerns. Moreover, the influence of the nonsymmetric electron distribution 
on the reconstructed result in MOT has not been extensively discussed and calls for more studies 
too. 

In this paper, the influence of large permanent dipoles on the MOT of nonsymmetric molecu- 
lar orbital is investigated. It is found that, the one-side-recollision condition in HHG can not be 
satisfied even with the few-cycle laser pulses, owing to the modification of the angle-dependent 
ionization rate induced by the Stark shift. We employ a tailored laser pulse to overcome this 
problem and the orbital reconstruction is carried out in both the length and velocity forms. The 
reconstructed results are compared to the ab initio orbital and the influence of the nonsymmetric 
electron distribution on the reconstructed result in MOT is discussed at the end. 

2. Theoretical model 

This section includes two independent parts: "Simulating the HHG" and "Reconstruction 
method". In the first part, we describe the modified Lewenstein model in the length form, which 



is employed in this work to simulate the HHG. In the second part, we introduce the orbital re- 
construction method in both the length and velocity forms. Note that in the second part the in- 
putting high-order harmonic spectra for the reconstruction process could be obtained from any 
ways: from calculations of different theories and forms or from experimental measurements. 

2.1. Simulating the HHG 

According to the Lewenstein model JT6| and using the approximate saddle points [22, 23 1, the 
time-dependent dipole moment is calculated by (atomic units are used throughout) 

D(f) = Y d aion{td)aprop(t d ,t)a rec [k(t d ,t)}+c.c., (1) 

td 

where a i0 „, a prop and a rec are the amplitudes of tunneling ionization, propagation after tunneling 
and recombination, t d is the approximate saddle time determined by solving 

~\p*(t d ,t)+A f (t d )] 2 = 0. (2) 

Af(t) is the vector potential corresponding to the electric field F(f), and p sf is the stationary 
momentum calculated by 

Pst(t d ,t) = / A f (t , )dt / . (3) 

t - td Jt d 

For atoms the amplitudes in Eq. (1) are given by 
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*ion{td) = j2T p exp[-^f^}, (4) 
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(5) 



a rec (k) = <y g \r\e ikr >. (6) 



Ip is the ionization energy and 1^ is the ground state in the form of hydrogenlike atom [ 16 1. 
k = p« + A/(f ) is the momentum of the return electron at the instant of recombination. 

For molecules, the ionization amplitude a,„„ is calculated with the MO-ADK theory ]24[ , 
which takes into account the structure features of molecular orbitals. *¥ g is the HOMO of the 



molecule obtained from Gaussian 03 ab initio code ]25) . For polar molecules whose HOMO 
has an unnegligible permanent dipole jj,/,, the ionization energy becomes time-dependent in 
response to the external electric field [20,2l][26): 



I p (t)=I p0 + n h -¥(t), (7) 
where I p o is the field free ionization energy. The permanent dipole /I/, is calculated by 

Hh = - fdrrp H (r). (8) 



p H (r) is the electron density of the HOMO calculated by p H (r) = ^(r)! 2 . 

For the Lewenstein model in the length form, a modification should be applied in Eq. (6) that 



the momentum k is replaced by the effective momentum [27 28 1 



(9) 



This modification implies that the recombination occurs within the potential well and the accel- 
eration effect of the binding potential is considered in this modified Lewenstein model p9][30| . 
In Eq. (9), the factor y(k) is given by 

1, k 2 >I p 



Y(k) = { , (10) 



sin 2 (ff). k 2 <I p 



This k-dependent factor avoids the unreasonable fact that all the return electrons have momenta 
larger than y/2I p . The introduction of y{k) also suppresses the strong background noises result- 
ing from the jump of k e /f over the gap from — ^j2~T p to \J2I P . As 7 = 1 for k 2 > I p , this term 
only affects the HHG of very low orders. 

Finally, the harmonic spectrum is obtained by Fourier transforming the time-dependent 
dipole acceleration a(f): 

Ejrt/v = Ja(t)exp(-iq(o L t)dt, (11) 
S g =| E xuv | 2 , (12) 

where a(f) = D(*), COl is the frequency of the driving pulse and q corresponds to the harmonic 
order. 

2.2. Reconstruction method 

As shown in Fig. 1, (x,y,z) is the molecular frame and (x',y',z') is the laboratory frame. The 
driving laser field is linearly polarized along the x' axis and propagates along the z' axis. To 
reconstruct the molecular orbital, one must calculate or measure the high-order harmonic spec- 
tra (including the intensity and phase) from the oriented target molecule at various orientation 
angles 0. In addition, one should also calculated or measure the high-order harmonic spectrum 
from the reference atom with the ionization energy similar to that of the target molecule. Then 
the recombination dipole moment d L and the dipole velocity d v for the molecular orbital in 
both x' and y' components are obtained by ]7) 

L/V 1 ET/U^^^M-i^Stirk^ 6 )} L/V 

d *>>^ = E^W^ d3) 

= 77(6,-) A%)j ^[^/(^^-^L^^-^W^/W, (14) 

where E, A, cp denote the complex electric field, amplitude and phase of the high-order har- 
monic radiation respectively. 77 is the scaling factor in molecules depending on the ionization 
angle 0, calculated by the MO-ADK theory, CO = qa>L is the harmonic frequency of the q-th 
order. The relation between CO and k is determined by the energy conversation CO = k 2 /2. To 

subtract the additional acquired Stark phase from E m ° l , the first-order Stark phase ^stark^ 03 ' ^) 
is introduced in Eqs. (13) and (14), which could be approximately evaluated with the analytical 



expression derived in 1 18 1. In the present paper, to focus our investigation on the influence of 
the permanent dipole in the aspect of modifying the ionization rate, the MOT is performed with 
the Stark phase perfectly subtracted in the simulation. In addition, the result of MOT with the 
Stark phase subtracted following JT8| as well as the influence of the permanent dipole in phase 
are briefly discussed at the end of Section 3. 

The dipole matrix elements of the "known" reference atom can be obtained by 

4 7 (k) = -id&„ f (k), (15) 
dj ef (k) = k^k); (16) 



where ^ re f(k) is the Fourier transform of the ground state of the reference atom *P re /(r). 




Fig. 1. Definition of the reference frames in the MOT scheme: (x,y,z) is the molecular frame 
and (x',y',z') is the laboratory frame. The driving laser field F is linearly polarized along 
the x' axis and propagates along the z' axis. 6 is defined as the angle between x and x' 
axes. 

By using the dipole moment or the dipole velocity, the reconstruction of the molecular orbital 
can be performed in the length form or in the velocity form. 
In the length form, the orbital can be reconstructed as: 

^(x,y) = -££4(a),0)exp{/[fe(a))cos(0)x + fe(a))sin(0)3;]}, (17) 
x e q 

*t*(x,y) = 1 £^^(Q),fl)exp{i[Jk(fl))cos(0>H-Jk(fl))sin(0)y]}, (18) 
y 6 q 

where 

4(co,6) =d%(<o,e)cos(6)-dy(a>,e)&m(e), (19) 
dy{(0,6) = rf£(a),e)sin(e)+<f£(a),0)cos(0). (20) 

W^(x,y) is the reconstructed orbital by using the dipole moments projected on the x axis of the 
molecular frame and *¥y {x,y) is the reconstructed orbital by using the dipole moments projected 
on the y axis. For symmetric molecules, the same orbital can in principle be reconstructed from 
both the x and y components of the dipole moments (d 1 ^ and dh ). However, due to the limited 
discrete sampling in Fourier space, they will most likely not give the same result [7|. Generally, 
the final reconstructed orbital is defined as 

^ L (x,y)= 1 -[^ L x (x,y)+^(x,y)} (21) 

to hopefully average out the distortions. 



Similarly, one can also reconstruct the orbital in the velocity form using the two projected 
components of the dipole velocity respectively: 




EE 



exp{i[k(a>) cos(9)x + k(co) sin(0)y]}, 



(22) 
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with 



4(a), 0) 



4£(ffl,0)cos(0)-<$(£O,0)sin(0), 
<^(a),0)sin(0)+^(a),0)cos(0). 



(24) 



(25) 



3. Result and discussion 

This section consists of two parts. In the first part, we analyze the HHG of a polar molecule 
to demonstrate the influence of the large permanent dipole on MOT and suggest a method to 
overcome this unfavorable influence. In the second part, we compare and discuss the recon- 
structed results in different forms and by using different projected components of the dipole 
matrix elements. 

3.1. Simulating and controlling the HHG 

In this work, we apply CO as the target molecule. The HOMO of CO has a large permanent 
dipole of 1 .72 a.u. pointing from C to O nucleus. The ionization energy of CO is 0.5 152 a.u., and 
the reference atom is Kr with the ionization energy of 0.5145 a.u. We apply a linearly polarized 
three-cycle sin 2 driving pulse F(t) = Fosin 2 (^)cos((»Lf + <f>o), where To is the optical cycle 
and <po = 1 .25 it. This form of laser pulse (with an opposite sign) has been used to accomplish 
the MOT of the nonsymmetric molecule HeH 2+ jl4J. The intensity and wavelength of the laser 
pulse we use are 1 x 10 14 W/cm 2 and 1.5 fxm respectively. Employing the 1.5 fim mid-IR laser 
source allows us to extend the harmonic spectrum considerably with a low laser intensity below 
the barrier suppression intensity | |T2"||3"T) , which ensures the validity of the used MO-ADK 
model and strong-field approximation (SFA) in our theoretical description. Another benefit of 
keeping the intensity low is that it will reduce the effect of the Stark shift. 

The electric field of the driving pulse is depicted by the dashed blue curve in Fig. 2(a). The 
HHG process is also illustrated by the black arrows in Fig. 2(a): under this three-cycle laser 
pulse, the HHG is dominantly contributed by the electrons ionized by the two peaks PI and 
P2. For the electrons ionized by PI, they first drift towards the positive direction and finally 
return to the core at times around Tl with negative momenta, emitting the high-order harmonic 
photons. On the contrary, the electrons ionized by P2 finally return to the core at times around 
T2 with positive momenta. 

When the effect of the permanent dipole of CO is ignored, the ionization energy is constant. 
The ionization probabilities in the +x' and -x' directions are dominantly determined by the 
waveform of the electric field, due to the exponential dependence of the tunneling ionization 
rate on the amplitude of electric field. Therefore, as the maximum amplitude of P2 is bigger than 
that of PI, electrons are more easily to be ionized by the P2 peak and the high-order harmonic 
radiation will be stronger at T2 than at Tl. This is confirmed by our calculation shown in the 
first row in Fig. (3). In Fig. 3(a), the red and green solid curves show the x' and y' components 
of the high-order harmonic spectrum from CO oriented at = 10° without taking into account 
the Stark shift. The harmonic spectrum from the reference atom is also shown in Fig. 3(a) by 




Fig. 2. (a) The electric fields of the driving, assistant and tailored pulses. The black ar- 
rows illustrate the two main contributing sets of trajectories in the HHG: the electron of 
PI— >T1 trajectory returns to the core from +x' position with a negative momentum, while 
the electron of P2— >T2 trajectory returns from the -x' side with positive momentum, (b) 
The time-dependent orbital energy Eg = — I p for = 10° and 170° respectively. 



the dashed curve. Figures 3(b) and (c) show the corresponding time-frequency analyses | 26||32) 
of the molecular HHG in the parallel (x') and perpendicular (y') directions respectively. It is 
found that the atomic spectrum exhibits double-plateau structure, while only the plateaus on 
the lower energy side are observed in the molecular spectra. This is because the plateaus on the 
higher energy side are so weak as to be below the background noise. The result implies that the 
plateaus on the lower energy side are at least 3-4 orders of magnitude higher than those on the 
higher energy side in both x' and y' components in the spectrum from CO. According to the 
time-frequency analyses, the plateaus on the lower energy side are mainly contributed by the Q2 
peaks, which are generated by the recollision of electrons with positive momenta at T2. While 
the plateaus on the higher energy side are contributed by the Ql peaks, which are generated 
by the recollision of electrons with negative momenta at Tl. From both the observations in the 
spectra (Fig. 3(a)) that the plateaus on the lower energy side are much higher than those on 
the higher energy side and in the time-frequency distributions (Figs. 3(b) and (c)) that Q2 are 
much stronger than Ql, it can be concluded that the high-order harmonics are generated by the 
recollision of electrons from only one side. 

However, when the stark shift is considered, the one-side-recollision condition is worse sat- 
isfied. As shown in Figs. 3(d), the plateaus on the higher energy side appear and are in the 
same level with the plateaus on the lower energy side in the harmonic spectra of CO in both 
components. The time-frequency analyses also show that the peaks Ql and Q2 have compara- 
ble intensities. It has been discussed that Ql is generated by the electrons with negative return 
momenta and Q2 is generated by the electrons with positive momenta. Therefore, the results 
indicate that the electrons with both negative and positive momenta have comparable contribu- 
tions to the HHG. 




Fig. 3. First row: (a) the x' and y' components of the high-order harmonic spectrum from 
the target molecule CO and (b)(c) the time-frequency distributions for each component. The 
HHG is calculated in response to the few-cycle pulse neglecting the effect of the permanent 
dipole. Second row: (d) the x' and y' components of the molecular harmonic spectrum and 
(e)(f) the time-frequency distributions for each component, after taking into account the 
Stark shift. The spectrum of the reference atom is also depicted by the dashed blue curve 
in panels (a) and (d). The values in the color bar correspond to lOZogio(Intensity). 



To explain why the Stark shift makes the one-side-recollision condition worse satisfied, the 
time-dependent orbital energy E g = —I p in response to the external electric field F{t) is plotted 
in Fig 2(b). The electron is more likely to be ionized with a high E g and is harder to be ionized 
when Eg is low. Therefore, the ionization rate of CO is decreased when /Z/, is parallel with the 
electric filed F and is increased in the antiparallel geometry. In the case of 10° orientation, as 
shown by the red solid curve in Fig. 2(b), the ionization by PI is enhanced by the Stark shift 
and the ionization by P2 is suppressed, leading to the comparable ionization probabilities from 
both peaks. And these comparable ionization probabilities finally lead to the comparable HHG 
contributions by electrons returning from opposite sides. 

As the MOT requires extracting structure information (the slices) of the target orbital from 
different directions, the one-side-recollision condition should be satisfied at various orientation 
angles. To evaluate the ratio between the probabilities that an electron returns to the core with 
negative and positive momenta at different angles, a semi -classical method introduced in |33) is 
used. First of all the electrons are tunneling ionized with the probability calculated by the MO- 
ADK theory. After ionization, the freed electron follows a classical trajectory, and we collect 
the electrons which finally return to the core. Every classical return receives a weight based 
on the ionization probability and a factor of T 3 , where T is the time the electron spends in the 
continuum from ionization to recombination and reflects the effect of wave-packet spreading. 
After the collection, the return probability of the electron with negative momentum P and 
the return probability of electron with positive momentum P + are obtained. The ratio of the 
probabilities P-/P+ for orientation angles from 0° to 180° is plotted in Fig. 4, where the blue 
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Fig. 4. The ratio between the probabilities that an electron returns to the core with negative 
and positive momentum for various orientation angles in the cases of: (1) using only the 
three-cycle driving pulse without taking into account the Stark shift (SS) (blue circles), (2) 
using only the three-cycle driving pulse with the SS (green diamonds), and (3) using the 
tailored pulse with the SS (red pentagrams). The inset shows the ratios on a logarithmic 
scale. 

circles and the green diamonds present the calculated results without and with the Stark effect 
taken into account respectively. Comparing the green diamonds to the blue circles, the Stark 
shift leads to more recollision with negative momentum especially for the orientation at small 
angles. The ratio of P / P + is bigger than 0. 1 for 9 < 60° and goes up to bigger than 1 for 9 = 
and 10°. 

To overcome this problem, one may first think of using an even lower intensity for the driving 
pulse to further reduce the effect of the Stark shift. However, too low intensity will bring the 
problem of a bad signal to noise ratio. Instead, in previous experimental and theoretical works 
higher intensities than 1 x 10 14 W/cm 2 were preferred, e.g. Ref. |12[|26| . Therefore, in this work 
we suggest a method to solve this problem by using a tailored laser pulse with the relatively 
high intensity. The tailored laser pulse is synthesized by adding a weak low-frequency laser 
field to the few-cycle laser field. The low-frequency field can be provided by different kinds 
of light sources such as the CO2 laser 1 34- 35] and the THz radiation |36][37}. The basic idea 
of adding this assistant field is to enlarge the difference between the amplitudes of PI and P2. 
In the present work, we apply the low-frequency pulse from CO2 laser with the wavelength of 
10.6 jJ.m. The intensity of this pulse is 1% of that of the few -cycle driving pulse and the electric 
field is F' — FQCO&(co' L t + (p' ) with <p' Q = l.ln. The low-frequency field and the tailored field 
are presented by the dash-dotted green and solid red curves in Fig. 2(a) respectively. 

The ratio P-/P+ under this tailored field is presented by the red pentagrams in Fig. 4. Com- 
paring to the green diamonds, the recollision with negative momentum is suppressed: the ra- 
tio is two orders of magnitude decreased and is below 0.1 for all the orientation angles (for 
180° < 9 < 360° the ratio at 9 equals that at 2k - 9). The result shows that the tailored laser 
pulse can efficiently restrict the recollision of the continuum electrons to only one side of the 
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Fig. 5. The calculated high-order harmonic spectra (the first column) and the corresponding 
time-frequency distributions for the x' and y' components of the molecular harmonic spec- 
tra (the second and third columns), by using the tailored pulse taking into account the Stark 
shift. The first, second and third rows correspond to 9 = 10°, 90° and 170° respectively. 
The values in the color bar correspond to 10/ogio(Intensity). 



core. 

The high-order harmonic spectra at three typical orientation angles 10°, 90° and 170° calcu- 
lated with the quantum model described in Section 2.1 are shown in the first column in Fig. 5. 
The second and third columns present the time-frequency analyses of the molecular harmonic 
spectra in the x' and y' components respectively. These results are consistent with those from 
the semi-classical method. The peaks Q2 are much stronger than Ql in all the time-frequency 
distributions. Regarding the molecular harmonic spectra, the plateaus on the lower energy side 
(predominantly contributed by electrons with positive momenta) are 3-4 orders higher than 
those on the higher energy side (contributed by electrons with negative momenta) at 9 = 10° as 
shown in Fig. 5(a). And the plateaus on the lower energy side are at least 4 orders of magnitude 
higher than those on the higher energy side for 9 = 90 and 170° as shown in Figs. 5(d) and 
(g). Therefore, the one-side-recollision condition is met and the HHG by this tailored pulse is 
satisfying for the MOT of the nonsymmetric target orbital. 

3.2. Tomographic reconstruction 

The MOT result in the length form is shown in Fig. (6). To reconstruct the orbital, we use 
the odd harmonics from the 23rd to the 79th orders (corresponds to 19-65 eV) generated from 
CO at 9 from 0° to 360° with the angular step A9 = 10°. Figures 6(a) and (b) display the 
reconstructed orbitals \P£ and \Pr using the dipole moments projected on the x and y axes in 
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Fig. 6. (a)(b) The orbitals reconstructed in the length forms by using the dipole moments 
projected on the x and y axis respectively. (c)(d) The intermediate functions f x (x,y) and 

fy{x,y). 

the molecular frame respectively. The ab initio orbital is also shown in Fig. 7(c). Compared 
with Fig. 7(c), Fig. 6(b) successfully reproduces the main structure of the nonsymmetric target 
orbital, while Fig. 6(a) fails to reproduce the structure. There is an artificial structure close to 
x=0 in Fig. 6(a) with a sharp jump from the negative value to the positive value. This is due to 
the nonsymmetric distribution of the orbital. 

To illustrate this, in Figs. 6(c) and (d) we present the intermediate functions f x and f y , which 
are defined as [6 1 

/^(^ J 3')=LL4( co ' ) ex p{ / ^( G) ) cos ( ) x+fc ( co ) sin ( ^^' 5 =x ^- (26) 

q 

The relation between X V L and / is 

Both f x (x,y) and f y (x,y) are symmetrically distributed along the y direction and are nonsym- 
metrically distributed along the x axis, corresponding to the nonsymmetric electric density dis- 
tribution of the target orbital along the permanent dipole. Therefore, one will have no problem 
in dividing f y by y to obtain because the nodal plane of f y locates right at y — 0. But the 
nodal plane is bent and not coincident with x = for f x . As a result, numerical problem arises 
when dividing f x by x with x— >0, which leads to the artificial structure and sharp jump in Fig. 
6(a). 




Fig. 7. (a)(b) The orbitals reconstructed in the velocity form by using the dipole velocities 
projected on the x and y axis respectively, (c) The ab initio orbital, (d) Slices of *Pf (LF), 
*Py (VF) and the ab initio orbital (AI) along the x axis. 



The reconstruction is also performed in the velocity form as shown in Figs. 7(a) and (b). The 
scan of 9 starts from 5° and increases in the step of 10° to avoid cos(0) = or sin(0) = 
in the denominators in Eqs. (22) and (23). *Fj and ^Jf are the reconstructed orbitals using the 
dipole velocities projected on the x and y axes respectively. Again, *Fj and give different 
results. WY successfully reproduces the main structure of the target orbital, while M'J is incor- 
rect. To further examine the influence of the nonsymmetric distribution of the orbital on the 
reconstructed result, we have also tried the MOT of another kind of orbital which is nonsym- 
metrically distributed along the y axis. On the contrary to the results shown in Figs. 6 and 7, 
and *F]f successfully reproduce the orbital this time, while "Ph and WY are incorrect (not 
shown). All the results indicate that, the MOT can only be accomplished by using the dipole 
matrix elements projected perpendicular to the permanent dipole in both the length and velocity 
forms for the nonsymmetric orbitals. 

Further comparing the reconstructed orbitals v Pr' and with the ab initio orbital, it is found 
that there is still detailed structure missed in the reconstructed results near the left core. This 
is due to the limited spectral range "detected" for the MOT fI7) . In Fig. 7(d), we compare 
the slices along the x axis of *¥y, and the ab initio orbital. The comparison shows that, 
the result obtained in the length form agrees with the ab initio orbital well for the positions and 
amplitudes of the left and middle peaks, but the position of the right peak deviates far away from 
the ab initio orbital. The result obtained in the velocity form matches the ab initio orbital better 
for the positions of all the three peaks, but the values at x=-3 and a.u. deviate a lot from those 
of the ab initio orbital. Both forms have limitation in reproducing the exact shape of the target 



orbital. Furthermore, the reconstruction of the orbital in the y direction is not satisfying. This 
error is also observed in the previous works especially for the a orbital. Besides the limited 
spectrum range, the imaging errors may also result from several other reasons, such as the 
discrete sampling in the frequency domain with the interval of 2o)l, and the SFA and single- 
active-electron approximation used in the theoretical model J7]- [T0}[T2) . In the MOT process, 
when the Stark phase is approximately subtracted by using the analytical expression derived 
in fTH) , the main structure of the target orbital is also successfully reproduced in both the length 
and velocity forms. The reconstructed orbitals have similar shapes to those shown in Figs. 6 and 
7, except that the whole structures are shifted towards the right in both forms. These shifts may 
therefore originate from the approximation in evaluating the Stark phase. On the other hand, 
this result implies that the influence of the permanent dipole in phase would lead to the shift of 
the reconstructed orbital, if the Stark phase is not exactly subtracted. 

4. Conclusion 

The influence of large permanent dipoles on the MOT is investigated in this paper. It is found 
that, owing to the modification of the angle-dependent ionization rate induced by the Stark 
shift of the orbital energy, the one-side-recollision condition can not be satisfied even with 
the few-cycle driving pulse. To overcome this problem, we employ a tailored laser pulse by 
adding a weak low-frequency pulse to the few-cycle pulse to control the HHG process and the 
requirement for the one-side-recollision is met. Then the reconstruction of the target orbital is 
performed in both the length and velocity forms. The results show that, the orbital structure can 
only be successfully reproduced by using the dipole matrix elements projected perpendicular 
to the permanent dipole in both forms. 
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